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Silicon materials play a key role in many technologically relevant fields, ranging from the elec¬ 
tronic to the photovoltaic industry. A systematic search for silicon allotropes was performed by 
employing a modified ab initio minima hopping crystal structure prediction method. The algo¬ 
rithm was optimized to specifically investigate the hitherto barely explored low-density regime of 
the silicon phase diagram by imitating the guest-host concept of clathrate compounds. In total 44 
metastable phases are presented, of which 11 exhibit direct or quasi-direct band-gaps in the range 
of 1.0-1.8 eV, close to the optimal Shockley-Queisser limit of ^ 1.4 eV, with a stronger overlap of 
the absorption spectra with the solar spectrum compared to conventional diamond silicon. Due to 
the structural resemblance to known clathrate compounds it is expected that the predicted phases 
can be synthesized. 


INTRODUCTION 

Silicon has become one of the technologically most rel¬ 
evant materials in the last century. It is the second most 
abundant element on earth and thus readily available, 
rendering it ideal for many industrial applications. Since 
elemental silicon is an intrinsic semiconductor that can 
be both p and n-doped it is well suited and thus widely 
used in the electronic industry, with applications in inte¬ 
grated circuits, photovoltaics and optoelectronic devices. 
Due to its dominance in electronics, it has been exten¬ 
sively studied both experimentally and theoretically such 
that its properties in the most stable form, the diamond 
structure {d-Si), are well understood. 

Although diamond silicon dominates the photovoltaic 
market due to its low cost, it is in fact not a good photo 
absorber. The band-gap of 1.12 eV in d-Si [1] is close 
to the maximal efficiency in the Shockley-Queisser limit 
of ^ 1.4 eV [2]. However, it is indirect with a large 
difference to the direct gap of 3.2 eV such that elec¬ 
tron excitations in the optical regime need to be phonon- 
mediated [3], preventing the use of d-Si in thin-film solar 
cells or other optical applications such as light-emitting 
diodes or high-performance transistors. Therefore, ef¬ 
forts have been made to design, discover and synthe¬ 
size direct band-gap materials with improved absorption 
properties. 

The phase diagram of silicon is complex with a se¬ 
quence of structures appearing in the high pressure 
regime [4, 5] such that several metastable phases of sili¬ 
con have been considered for use in photovoltaic applica¬ 
tions [6-9]. Allotropes such as the BC8, R8 or londs- 
daleite silicon however turn out to be unsuited since 
they are either semi-metallic (BC8), have a too small 
(0.24 eV for R8) or indirect band-gaps (^1 eV for londs- 
daleite) [8, 10]. Recently, novel silicon allotropes have 


been proposed and studied theoretically through struc¬ 
tural search methods, such as the Si20T [11] and P- 

I [12] phases, which exhibit promising electronic prop¬ 
erties with either direct or quasi-direct band-gaps in the 
desired energy range. Carefully tailored super-lattices of 
interlaced silicon and germanium layers were also pro¬ 
posed by D’Avezac et al. [13] with excellent absorption, 
but techniques for the required atomic scale layer-by- 
layer assembly on an industrial scale are yet to be de¬ 
veloped. 

On the other hand, the low-density regime of the sil¬ 
icon phase diagram remains relatively poorly explored, 
although semi-conducting group IV clathrates [14-16] 
have recently drawn attention as a possible material 
for various energy applications. Since the structure of 
clathrates is characterized by face-sharing cage-like units 
of host atoms enclosing guests, its properties can be 
tuned through tailored substitution of the host or guest 
elements. Superconductivity was observed in clathrate 
compounds [17, 18] with transition temperatures as high 
as 8 K in the barium-encapsulated type I silicon clathrate 
BagSUe synthesized at high pressure [19]. Further¬ 
more, the properties induced by the vibrational “rat¬ 
tling” modes of the guest atoms in the large host cages 
leads to reduced thermal conductivity at high electronic 
conductivity (the phonon glass-electron crystal concept 
introduced by Slack [20]), such that clathrate have been 
studied in the context of thermoelectric applications [21- 
26]. More recently, these materials have also been investi¬ 
gated for their potential use in lithium ion batteries [27- 
29] and as hydrogen storage materials [30] due to the 
large voids that could be charged/discharged with guest 
atoms or molecules. 

Guest free group IV element clathrates of type I and 

II have been touted as promising candidates for photo¬ 
voltaic applications due to the direct band-gap character 
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with gap energies close to the optimal value [31-33]. The 
type II silicon clathrate has a wide direct band-gap of 
1.9 eV, which is however not allowed for electronic dipole 
transitions [25]. Nevertheless, Baranowski et al investi¬ 
gated the possibility of tuning the electronic properties 
of the type II clathrate through variation of the chem¬ 
ical composition in Si-Ge alloys, resulting in band-gaps 
between 0.8-1.8 eV [31]. Furthermore, novel experimen¬ 
tal techniques were developed to efficiently prepare large 
samples of phase-pure type II silicon clathrates, an essen¬ 
tial requirement for future large-scale industrial develop¬ 
ments [33]. 

Although only group IV clathrates of type I, II, III, 
VIII and IX were observed experimentally, theoretical 
studies indicate that there is a large structural diver¬ 
sity of low-density silicon phases [34]. Structural prin¬ 
ciples of various polytypes and intergrowth frameworks 
of clathrates were systematically investigated from basic 
and extended polyhedral building blocks by Karttunen et 
al. [35], while solutions to the Kelvin problem, the task 
of partitioning three-dimensional space into cells of equal 
volume with minimal area, were used to predict novel 
low-density silicon phases [36]. More recently, global 
search methods were used to predict silicon allotropes 
with possible applications in optoelectronics. Wang et 
al. [37] proposed 6 different metastable silicon allotropes 
exhibiting channel-like structures with direct or quasi- 
direct band gaps predicted with particle swarm optimiza¬ 
tion, and a novel cage-like structure was discovered by 
Nguyen et al. [38] by a genetic algorithm. 

Experimentally, all routes to synthesize guest-free low 
density silicon phases involve, in some way or another, 
precursor materials containing silicon and guest elements. 
Thermal decomposition of Zintl monosilicide phases is 
one possible route to obtain silicon clathrates [39, 40]. 
In this process, the cations sublimate, causing a charge 
imbalance, and leading to the four-fold rearrangement of 
the silicon atoms around the remaining template cations. 
Other methods include chemical oxidation, spark plasma 
treatment of Na 4 Si 4 [41] or directly from elemental Si 
and Na at high pressure [42, 43]. The guest elements are 
then thermally removed, a process referred to as “thermal 
degassing”. In fact, with this procedure it is possible to 
synthesize metastable silicon polymorphs with energies 
considerably beyond the ground state [43]. 

In this work we conduct a computational search for 
novel low-density silicon allotropes with potential use in 
photovoltaic applications. Using first principle calcula¬ 
tions, we discover a plethora of hitherto unknown low- 
density low-energy silicon phases consisting of interlinked 
cage-like polyhedra and channel-like structures. 


METHOD 

Inspired by the existence of group IV clathrates, the 
minima hopping crystal structure prediction method 
(MHM) [44, 45] was adapted to investigate low-density 
silicon allotropes. The MHM was initially designed to 
thoroughly scan the low-lying energy landscape of any 
material and to identify stable and metastable phases 
by performing consecutive short molecular dynamics es¬ 
cape steps followed by local geometry relaxations taking 
into account both atomic and cell variables. The poten¬ 
tial energy surface is mapped out efficiently by aligning 
the initial molecular dynamics velocities approximately 
along soft mode directions [46], thus exploiting the Bell- 
Evans-Polanyi [47, 48] principle to steer the search to¬ 
wards low energy structures. The predictive power of 
this approach has been demonstrated in a wide range of 
applications [47, 49-52]. 

The structural search was performed using the Lenosky 
tight-binding scheme to model the Si-Si interaction [53], 
which has shown to give sufficiently accurate results for 
solids at moderate computational cost [54, 55]. To im¬ 
itate encapsulated guest elements, fictitious Lennard- 
Jones (LJ) spheres were placed in the simulation cell as 
inert template guest atoms. The LJ potential was shifted 
and truncated such that it decreases continuously to zero 
at a cutoff radius Vc and no discontinuities arise in the 
forces [56]. To avoid the formation of LJ dimers the in¬ 
teraction between the LJ particles was further truncated 
at the minimum and shifted towards zero, resulting in a 
purely repulsive potential. The parameters for the LJ- 
silicon interaction were carefully tuned to already known 
silicon clathrate structures. In analogy to experimentally 
used inert template elements such as Na, the interaction 
between the LJ-sphere and the silicon framework was de¬ 
signed to be minimal, with a strong repulsive potential 
and a weak attraction. The parameters were thus set to 
a = 3.00 A, e = 0.02 eV, Tc = 7a. 

Density functional theory calculations were carried out 
to refine the structures with the projector augmented 
wave formalism as implemented in the VASP package [57] 
using the Perdew-Burke-Erzernhof (PBE) [58] exchange- 
correlation functional. A plane-wave cutoff energy of 
500 eV was used together with a sufficiently dense k-point 
mesh, resulting in total energies converged to less than 
I meV/atom. Both atomic and cell variables were relaxed 
simultaneously until the forces on the atoms were less 
than 3 meV/A and the stresses were less than O.I eV/A^. 

Electronic properties were investigated by calculat¬ 
ing the Kohn-Sham electronic band structure within the 
PBE approximation to the exchange-correlation poten¬ 
tial. Since conventional DET systematically underesti¬ 
mates the band-gap, the hybrid HSE06 functional [59-61] 
was used on a set of selected structures including d-Si and 
the type II clathrate. The band-gap error within PBE 
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was found to lie consistently between 0.4 and 0.65 eV. 
Therefore, a scissor operator was subsequently used to 
obtain a corrected gap by shifting the unoccupied bands 
by an averaged value of 0.55 eV [62]. 

Absorption spectra were obtained within time- 
dependent density-functional theory [63, 64] using the 
code YAMBO [65] starting from PBE ground-states ob¬ 
tained with ABINIT [66]. We used the semi-empirical 
ajc^ long-range exchange-correlation kernel [67], that 
was shown to give very good results for systems with 
delocalized Wannier-Mott excitons (such as the ones we 
expect in Si-based materials) [67, 68]. We used the (con¬ 
stant) value of a = —0.2, which is the one determined for 
bulk silicon. For the optical spectra we used randomly 
shifted /c-point grids with the number of /c-points in di¬ 
rection a equal to 71 x jbc^j, where b is the reciprocal 
lattice vector in direction a. The number of bands in¬ 
cluded in the calculation was equal to 5 times the num¬ 
ber of atoms, i.e., we used 3 unoccupied bands per Si 
atom. This amounts toal2xl2xl2 /c-point mesh and 
6 unoccupied bands for d-Si. 

The phonon calculations were carried out with 
the frozen phonon approach as implemented in the 
PHONOPY package [69]. Super-cells of dimension 4 x 4 x 
4 were used for cells with less than 12 atoms, 3 x 3 x 3 for 
cells with between 13 and 24 atoms, and 2x2x2 super¬ 
cells for larger unit cells, which were found to give well 
converged results. The X-ray diffraction pattern were 
simulated with MERCURY [70] with the Cu Ka radia¬ 
tion of wavelength 1.54056 A. 

RESULTS AND DISCUSSION 

The MHM simulations at the tight-binding level were 
performed for different stoichiometries Si^^LJ^, where x 
and y were chosen from two different approaches. First, 
a systematic set of stoichiometries for moderate system 
sizes were explored with all combinations of x and y in 
the range oi x = 16,18,..., 34 and ^ = 2,4, respectively. 
The initial structures were generated randomly, without 
any symmetry constraints. Second, a set of known struc¬ 
tures from literature was used to generate initial struc¬ 
tures for the MHM simulations. The energetically lowest 
energy clathrates of type I and II and structures of the 
Kelvin-type from Ref. [36] were used. For each structure, 
the center of the polyhedra was located and LJ spheres 
were placed accordingly. The stoichiometries [x\y] then 
resulted in [4618] from Type I, [3416] from Type II, [1212] 
from Si{^, [4618] from Si{^, [82|14] from Si{^^, [4117] from 
Si^y, [4017] from Si^, [178130] from Si^^, [68|12] from 
Si^^^ [58110] from [4618] from Si{^ and [232140] 

from Si^. 

Each MHM run resulted in several thousand struc¬ 
tures, of which all with energies (calculated in tight- 
binding) 400 meV/atom higher than d-Si were discarded. 


Subsequently, the LJ spheres were removed and the re¬ 
maining structures were carefully filtered for duplicates 
and analyzed with respect to symmetry. Due to the 
large size of the system under investigation, the ma¬ 
jority of the phases encountered during the structural 
search were amorphous or some defect structures of a 
crystalline phase. In fact, less than 5% of the phases 
encountered during a search had at least one symmetry 
element. Therefore, all structures that had no symme¬ 
try at all were discarded, unless they were very low in 
energy. Finally, the energetically most promising crys¬ 
talline phases were refined at the DFT level. Both the 
atomic and cell variables were relaxed to the nearest lo¬ 
cal minimum without constraining the initial symmetry 
to avoid relaxation to a first or higher-order saddle point. 

After the initial filtering, roughly 150 structures were 
considered for further analysis. Among those, several 
structures were rediscovered that had already been dis¬ 
cussed in literature, such as the bonded chiral framework 
structure (CFS), first reported by Pickard et a/., which 
is energetically degenerate to the type H clathrate. Al¬ 
though the structure does not contain any cages, the vol¬ 
ume per atom is considerably higher than d-Si with a 
value of 22.17 due to the channel-like voids. Similarly, 
a structure with P4:2/mnm symmetry was found that was 
reported earlier by Nguyen et al. [38] with an energy of 
114 eV/atom above d-Si and a volume of 23.31 A^. The 
H-clathrate, initially described by Udachin et al [71], was 
also found with an energy of 83 eV/atom above d-Si and 
a volume of 24.09 A^. 

FIG. 1. Hybrid channels embedded in diamond structure with 
atoms denoted by yellow (light gray) spheres, (a) Structural 
motif representing dimers (red or dark gray spheres) within a 
8-ring channels (green or medium gray) and (b) 4-rings (red or 
dark gray spheres) within 10-ring channels (green or medium 
gray). 


(a) (b) 

As expected, all structures are sp^-bonded and have 
lower densities than d-Si due to the fictitious LJ spheres. 
In general, the low density in all allotropes is achieved 
by one of two different structural motifs. The first type 
consist of wide 1-dimensional channels, which either run 
in one, two or three dimensions, can be straight, or even 
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run in zig-zag lines through the material. Some represen¬ 
tatives of this class have been extensively studied in the 
past for various sp^ materials. It has been shown that 
low-energy low-density allotropes with 5- to 8 -membered 
rings can be obtained systematically from different stack¬ 
ing of cubic diamond layers, and structures including 
combinations of 4-, 6 - and 8 -rings have been reported 
in various publications [ 12 , 34, 52, 72-74]. Energetically, 
these structures are close to d-Si, and the channels can 
exist as defect wires or planes embedded in the cubic 
or hexagonal diamond structure [73]. The recently syn¬ 
thesized open framework phase of Kim et al [43], which 
contains both 5- and 8 -membered rings, can be classified 
to belong to this type and was also discovered during 
our search. Similar to clathrate structures, the channels 
are large enough to host guest atoms [42]. Since these 
phases have been already widely discussed in literature, 
all structures were put into the supplementary materials. 
We however also discovered two other channel-like mo¬ 
tifs which frequently appeared in low-energy allotropes 
during our search, namely the hybrid channels enclosing 
dimers within 8 -rings or enclosing 4-rings within 10-rings, 
as shown in Figure 1. The latter can be interpreted as 
chains of [4^5^6jjj] cages within a diamond framework 
(see below). 

In contrast to the first class of structures, the second 
contains fully 3-dimensional cages of different types and 
sizes as shown in Figure 2 . Similar to the known clathrate 
structures, the basic building blocks are polyhedra with a 
number of silicon atoms ranging from 14 to 40, consisting 
of faces with 4, 5, 6 , 7 and 8 -membered rings. Each poly¬ 
hedron can be described by a set of indexes with corre¬ 
sponding superscripts in square brackets, [A^B^C^ ... ], 
where ^4, 5, C, ... correspond to the number of ver¬ 
texes of each face, and a, /3, 7 , ... correspond to their 
multiplicity. Thus, the type I clathrate would be de¬ 
scribed as [ 5 ^^] 2 [ 5 ^^ 6 ^] 6 , indicating that the structure is 
composed of 2 pentagonal dodecahedra and 6 hexagonal 
truncated trapezohedra. Overall, cage-like voids lead to 
significantly lower densities than channels in the mate¬ 
rial. 

Figure 3 and Table I show a selection of 44 low-energy 
allotropes consisting of cage-like voids which were char¬ 
acterized in detail. Most structures are made from inter¬ 
linked polyhedral building blocks of 41 different kinds, 
which can be categorized into polyhedra with planar 
or nearly planar faces (first two panels in Figure 2 ) or 
with at least one curved face (three last panels in Fig¬ 
ure 2). The stability of the individual cages was evalu¬ 
ated in vacuum with the BigDFT wavelet code [75] to¬ 
gether with Hartwigsen-Goedecker-Hutter pseudopoten¬ 
tials [76]. The calculations were performed imposing zero 
boundary conditions on the wave-functions and the trun¬ 
cated bonds were passivated with hydrogen atoms. Cages 
containing 4-rings were found to be energetically unfa¬ 
vorable due to the strong strain on the bond angles, the 


cohesive energy showing a strong correlation with the 
number of 4-rings in the cages, indicated by a high cor¬ 
relation coefficient of 0.89. On the other hand, the pen¬ 
tagonal dodecahedra [5^^] is the least strained and thus 
energetically favorable. To some extent, this trend car¬ 
ries on to the solid, where the content of 4-rings increases 
with the energy of the phases with a correlation coeffi¬ 
cient of 0.57. On the other hand, [5^^] appears as of¬ 
ten as 15 times among 44 phases. Two phases, namely 
C30 {Pna2‘l^) and C42 /-4^^, could not be clearly de¬ 
composed into polyhedra, although they have an open 
framework structure of sp^ bonded silicon atoms. 

Phonon calculations were carried out to assess the dy¬ 
namical lattice stability for every phase. No imaginary 
phonon modes were found, thus confirming their struc¬ 
tural stability (see supplemental materials). From the 
energetic point of view, all structures are less favorable 
than d-Si or the most stable silicon clathrate of type II 
and are thus metastable. The relative energies with re¬ 
spect to d-Si are plotted as a function of the atomic vol¬ 
ume for all 44 phases in Figure 4. However, at least the 
20 lowest energy structures ranging from COl to C 20 are 
below 150 meV/atom in energy relative to d-Si and are 
thus likely to be thermodynamically accessible in exper¬ 
iments. Furthermore, 8 of those structures have energies 
below or very close to the d-Si line at their specific equi¬ 
librium volumes, namely COl, C03, C04, 006, C07, C08, 
CIO and C13, thus representing the most promising set 
of candidates for future experimental synthesis. Further¬ 
more, the bulk moduli Bq were computed through a fit 
to the Murnaghan equation of states (see table I) and are 
consistently lower than the value for d-Si of 87.9 GPa. 

Clathrate compounds are characterized by a frame¬ 
work of an open host structure with cavities occupied 
by guest molecules or atoms, and naturally the size of 
the cavities determines the maximum size of the guest 
atoms hosted by the crystal lattice. Although Na is 
the most prominent guest atom in silicon clathrates, 
clathrate materials can host various guest atoms which 
interact weakly with the host atoms based on the Zintl 
concept and profoundly change the electronic proper¬ 
ties [14, 19]. To asses the size of guest molecules that 
can be hosted by the cavities, we estimated the radii of 
the largest possible spheres placed at the center of each 
polyhedra without distorting the host structure. All sil¬ 
icon atoms were assumed to be rigid spheres with a co¬ 
valent radius of 1.11 A. Since the average bond-length in 
the solid is slightly larger than twice the covalent ra¬ 
dius (^2.38 A) and all atoms are fully coordinated with 
no bonds pointing towards the center of the voids, we 
can safely assume that this value is an upper limit. The 
smallest cavity of type [4^5^] can thus host guests with 
a radius of merely 1.3 A, whereas the largest could en¬ 
close molecules with radii of 3.3 A corresponding to the 
[ 4258612 ] cage. The average guest radius for all structures 
is listed in table I is 2.0 A. Nevertheless, the estimated 
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FIG. 2. Polyhedral cages which constitute the building blocks for the crystal lattice of the novel silicon allotropes. 


[4356] [4258] [435^6^] [4‘‘5‘‘63] [435^6^] [5“] [415*06^] [4®6®] 



[5“62] [5“6®] [4‘'5''6«] [S'^e^] [4^5106®] [5'^6®] [425861^] 



[42586^] [4^6®] [42586^] [4^5'‘6'‘] [4V8^] [485^8^] [4^586^ [48586%] 



[48586878] [586% [48586% [48586^1] [48586%] [48586^1,] [485868] [485^68] 



[485868] [4‘‘5''68] [4^5868] [485868] [4^5V„] [4'858068] [58858] [58858] [5885^78] 
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spherical radii for the guest atoms should be considered 
with care and should only be taken as a guide for the 
maximal size of inert guest molecules. 

Most structures are semi-conducting with gaps in the 
range between about 0.6 and 2.0 eV as shown in table I. 
Several structures exhibit a direct or quasi-direct band 
gap, namely COl, COS, C06, COT, COS, C12, C28, C32, 
CSS, C35 and C36 (see table I and 103, Ill, 112, 124, 142, 
146 in the supplemental material). The most promising 
candidates for photovoltaic applications are thus COS, 
COT, COS, C32 and C36, which have gaps close to the 
optimal value for the Shockley-Queisser limit of ~ 1.4 eV. 

In Figure 5 the calculated absorption spectra of some of 
the most promising allotropes are shown. The calculated 
absorption spectrum of the cubic diamond phase is also 
shown for comparison, together with the reference air 
mass (AM) 1.5 solar spectral irradiance [TT]. We see that, 
contrary to d-Si, most of these new allotropes absorb well 
in the visible, in a region perfectly compatible with the 
solar spectrum. This is true even if only a few structures 
have a direct or quasi-direct band gap. The absorption 
in this region is still below the ones of the Si20T [11] 
and P-1 [12] phases proposed before, but still sufficient 
to allow for thin-film Si photovoltaic cells. 


CONCLUSION 

In conclusion, a systematic search was performed using 
a modified version of the minima hopping method, tai¬ 
lored to imitate the synthesis process of silicon clathrates, 
to identify low density silicon allotropes. The efficiency 
of our method was confirmed by the discovery of many 
phases reported in literature and a plethora of hitherto 
unknown phases. The two main structural motifs lead¬ 
ing to low density are the formation of channels or cages. 
Due to the large sizes of the voids both in the chan¬ 
nels and cages they can readily host guest atoms which 
can be used to tune the material properties. Further¬ 
more, such guest-host materials are likely to be stabilized 
by pressure as demonstrated by Kurakevych et al. [42]. 
In view of the potential use in optical applications the 
electronic properties were investigated by computing the 
band structure and the absorption spectra. Several low- 
density allotropes show a significantly improved overlap 
with the solar spectrum as compared to the conventional 
diamond silicon and are thus promising candidates for 
use in thin-film photovoltaic applications. 

This work was done within the NCCR Marvel project. 
MA gratefully acknowledges support from the Novar¬ 
tis Universitat Basel Excellence Scholarship for Life Sci¬ 
ences. We acknowledge the computational resources 
provided by the Swiss National Supercomputing Center 
(CSCS) in Lugano (project s499) and GENCI (project 
x201409601T) in Erance. 


* Stefan.Goedecker@unibas.ch 

[1] W. Bludau, A. Onton, and W. Heinke, Journal of Ap¬ 
plied Physics 45, 1846 (19T4). 

[2] W. Shockley and H. J. Queisser, Journal of Applied 
Physics 32, 510 (1961). 

[3] P. Lautenschlager, M. Garriga, L. Vina, and M. Gar- 
dona, Phys. Rev. B 36, 4821 (198T). 

[4] S. Duclos, Y. Vohra, and A. Ruoff, Phys. Rev. B 41, 
12021 (1990). 

[5] M. I. McMahon and R. J. Nelmes, Phys. Rev. B 47, 833T 
(1993). 

[6] B. D. Malone and M. L. Gohen, Phys. Rev. B 85, 024116 

( 2012 ). 

[T] B. D. Malone, S. G. Louie, and M. L. Gohen, Phys. Rev. 
B 81, 115201 (2010). 

[8] B. D. Malone, J. D. San, and M. L. Gohen, Phys. Rev. 
B 78, 161202 (2008). 

[9] B. D. Malone, J. D. San, and M. L. Gohen, Phys. Rev. 
B 78, 035210 (2008). 

[10] J. M. Besson, E. H. Mokhtari, J. Gonzalez, and G. Weill, 
Phys. Rev. Lett. 59, 473 (1987). 

[11] H. J. Xiang, B. Huang, E. Kan, S.-H. Wei, and X. G. 
Gong, Phys. Rev. Lett. 110, 118702 (2013). 

[12] S. Botti, J. A. Flores-Livas, M. Amsler, S. Goedecker, 
and M. A. L. Marques, Phys. Rev. B 86 , 121204 (2012). 

[13] M. d’Avezac, J.-W. Luo, T. Ghanier, and A. Zunger, 
Physical Review Letters 108, 027401 (2012). 

[14] D. M. Rowe, Thermoelectrics Handbook: Macro to Nano, 
1st ed. (GRG Press, 2005). 

[15] M. Beekman and G. S. Nolas, J. Mater. Ghem. 18, 842 
(2008). 

[16] K. A. Kovnir and A. V. Shevelkov, Russian Ghemical 
Reviews 73, 923 (2004). 

[17] H. Kawaji, H.-o. Horie, S. Yamanaka, and M. Ishikawa, 
Phys. Rev. Lett. 74, 1427 (1995). 

[18] Y. Li, J. Garcia, N. Ghen, L. Liu, F. Li, Y. Wei, S. Bi, 
G. Gao, and Z. S. Feng, Journal of Applied Physics 113, 
203908 (2013). 

[19] S. Yamanaka, E. Enishi, H. Fukuoka, and M. Yasukawa, 
Inorg. Ghem. 39, 56 (2000). 

[20] G. S. Nolas, J. L. Gohn, G. A. Slack, and S. B. Schujman, 
Applied Physics Letters 73, 178 (1998). 

[21] G. Nolas and G. Slack, American Scientist 89, 136 (2001). 

[22] P. Rogl (IEEE, 2005) pp. 440-445. 

[23] A. San-Miguel and P. Toulemonde, High Pressure Re¬ 
search 25, 159 (2005). 

[24] A. Saramat, G. Svensson, A. E. G. Palmqvist, G. Stiewe, 
E. Mueller, D. Platzek, S. G. K. Williams, D. M. Rowe, 
J. D. Bryan, and G. D. Stucky, Journal of Applied 
Physics 99, 023708 (2006). 

[25] D. Gonnetable, Phys. Rev. B 75, 125202 (2007). 

[26] M. Ghristensen, S. Johnsen, and B. B. Iversen, Dalton 
Trans. 39, 978 (2010). 

[27] J. Yang and J. S. Tse, J. Mater. Ghem. A 1, 7782 (2013). 

[28] T. Langer, S. Dupke, H. Trill, S. Passerini, H. Eckert, 

R. Pttgen, and M. Winter, J. Electrochem. Soc. 159, 
A1318 (2012). 

[29] N. A. Wagner, R. Raghavan, R. Zhao, Q. Wei, X. Peng, 
and G. K. Ghan, GHEMELEGTROGHEM 1, 347 (2014). 

[30] N. Okamoto, D. Neiner, Q. Ramasse, G. Gondron, 

S. Kauzlarich, and N. Browning, Microscopy and Mi- 



7 


FIG. 3. Low-density silicon allotropes formed of cage-like building blocks. The colors of the polyhedra correspond to the ones 
in Figure 2. 



F 222128 | 


Fmmm"® 


Amm2^^^ 




Cmcm^^ 




Pmma'^^ 


Ama2^ 


[ 5 ’ 7 » 8 M 


[4‘5‘6‘8n 


[4*5«6’8il 


(4‘6''83] 


[324I52668I] 


Pmn2F°ll 


Pmn2, 


Pmmm'^° 


Pbcn^^ 


Pmmmi^^ 


Cccm^ 


[4^5‘'6*8^] 


(4^5V8^] 


iAmm2^° 


Ima2^^ 


Pnma^^ 


t4^5‘'6='7^1 


[435^6^7^] 


(5^6''8'] 


[s'e-'sM 


Pnma- 


Cmmm^^ 


Immm' 


[4^5®72] 


[4^6S7^] 


[425*6^72] 


Immm' 











croanalysis 13, 1176 (2007). 

[31] L. L. Baranowski, L. Krishna, A. D. Martinez, T. Ra- 
harjo, V. Stevanovic, A. C. Tamboli, and E. S. Toberer, 

J. Mater. Chem. C 2, 3231 (2014). 

[32] A. Martinez, L. Krishna, L. Baranowski, M. Lusk, E. To¬ 
berer, and A. Tamboli, IEEE Journal of Photovoltaics 
3, 1305 (2013). 

[33] L. Krishna, L. L. Baranowski, A. D. Martinez, C. A. 
Koh, P. C. Taylor, A. C. Tamboli, and E. S. Toberer, 
CrystEngComm 16, 3940 (2014). 

[34] M. A. Zwijnenburg, K. E. Jelfs, and S. T. Bromley, Phys. 
Chem. Chem. Phys. 12, 8505 (2010). 

[35] A. J. Karttunen, T. E. Eassler, M. Linnolahti, and T. A. 
Pakkanen, Inorg. Chem. 50, 1733 (2011). 

[36] H.-Y. Zhao, J. Wang, Q.-M. Ma, and Y. Liu, Phys. 
Chem. Chem. Phys. 15, 17619 (2013). 

[37] Q. Wang, B. Xu, J. Sun, H. Liu, Z. Zhao, D. Yu, C. Pan, 
and J. He, J. Am. Chem. Soc. 136, 9826 (2014). 

[38] M. C. Nguyen, X. Zhao, C.-Z. Wang, and K.-M. Ho, 
Phys. Rev. B 89, 184112 (2014). 

[39] J. Gryko, P. E. McMillan, R. E. Marzke, G. K. Ra- 
machandran, D. Patton, S. K. Deb, and O. E. Sankey, 
Phys. Rev. B 62, R7707 (2000). 

[40] M. Beekman, C. P. Sebastian, Y. Grin, and G. S. Nolas, 
Journal of Elec Materi 38, 1136 (2009). 

[41] M. Beekman, M. Baitinger, H. Borrmann, W. Schnelle, 

K. Meier, G. S. Nolas, and Y. Grin, J. Am. Chem. Soc. 
131, 9642 (2009). 

[42] 0.0. Kurakevych, T. A. Strobel, D. Y. Kim, T. Mura- 
matsu, and V. V. Struzhkin, Crystal Growth Sz Design 
13, 303 (2013). 

[43] D. Y. Kim, S. Stefanoski, 0.0. Kurakevych, and T. A. 
Strobel, Nat Mater 14, 169173 (2015). 

[44] S. Goedecker, J. Chem. Phys. 120, 9911 (2004). 

[45] M. Amsler and S. Goedecker, J. Chem. Phys. 133, 224104 

( 2010 ). 

[46] M. Sicher, S. Mohr, and S. Goedecker, The Journal of 
Chemical Physics 134, 044106 (2011). 

[47] S. Roy, S. Goedecker, M. J. Eield, and E. Penev, J. Phys. 
Chem. B 113, 7315 (2009). 

[48] E. Jensen, Introduction to Computational Chemistry: 
Second Edition (JW, 2011). 

[49] W. Hellmann, R. G. Hennig, S. Goedecker, C. J. Umri- 
gar, B. Delley, and T. Lenosky, Phys. Rev. B 75, 085411 
(2007). 

[50] K. Bao, S. Goedecker, K. Koga, E. Lancon, and 
A. Neelov, Phys. Rev. B 79, 041405 (2009). 

[51] A. Willand, M. Gramzow, S. A. Ghasemi, L. Genovese, 
T. Deutsch, K. Reuter, and S. Goedecker, Phys. Rev. B 
81, 201405 (2010). 

[52] M. Amsler, J. A. Elores-Livas, L. Lehtovaara, E. Bal- 
ima, S. A. Ghasemi, D. Machon, S. Pailhes, A. Willand, 
D. Caliste, S. Botti, A. San Miguel, S. Goedecker, and 
M. A. L. Marques, Phys. Rev. Lett. 108, 065501 (2012). 

[53] T. J. Lenosky, J. D. Kress, 1. Kwon, A. E. Voter, B. Ed¬ 
wards, D. E. Richards, S. Yang, and J. B. Adams, Phys. 
Rev. B 55, 1528 (1997). 

[54] S. A. Ghasemi, M. Amsler, R. G. Hennig, S. Roy, 

S. Goedecker, T. J. Lenosky, C. J. Umrigar, L. Genovese, 

T. Morishita, and K. Nishio, Phys. Rev. B 81, 214107 


( 2010 ). 

[55] S. A. Ghasemi, T. J. Lenosky, M. Amsler, A. Sadeghi, 

L. Genovese, and S. Goedecker, Phys. Rev. B 90, 054117 
(2014). 

[56] S. D. Stoddard and J. Eord, Phys. Rev. A 8, 1504 (1973). 

[57] G. Kresse and E. J., Comput. Mat. Sci. 6, 15 (1996). 

[58] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

[59] J. Heyd, G. E. Scuseria, and M. Ernzerhof, J. Chem. 
Phys. 118, 8207 (2003). 

[60] J. Paier, M. Marsman, K. Hummer, G. Kresse, 1. C. Ger¬ 
ber, and J. G. Angyan, J. Chem. Phys. 125, 249901 
(2006). 

[61] J. Heyd, G. E. Scuseria, and M. Ernzerhof, J. Chem. 
Phys. 124, 219906 (2006). 

[62] Z. H. Levine and D. C. Allan, Phys. Rev. Lett. 63, 1719 
(1989). 

[63] M. Marques and E. Gross, Annual Review of Physical 
Chemistry 55, 427 (2004). 

[64] S. Botti, A. Schindlmayr, R. D. Sole, and L. Reining, 
Reports on Progress in Physics 70, 357 (2007). 

[65] A. Marini, C. Hogan, M. Grning, and D. Varsano, Com¬ 
puter Physics Communications 180, 1392 (2009). 

[66] X. Gonze, B. Amadou, P.-M. Anglade, J.-M. Beuken, 
E. Bottin, P. Boulanger, E. Bruneval, D. Caliste, R. Cara¬ 
cas, M. Ct, T. Deutsch, L. Genovese, P. Ghosez, M. Gi- 
antomassi, S. Goedecker, D. Hamann, P. Hermet, E. Jol¬ 
iet, G. Jomard, S. Leroux, M. Mancini, S. Mazevet, 

M. Oliveira, G. Onida, Y. Pouillon, T. Rangel, G.-M. 
Rignanese, D. Sangalli, R. Shaltaf, M. Torrent, M. Ver- 
straete, G. Zerah, and J. Zwanziger, Computer Physics 
Communications 180, 2582 (2009). 

[67] L. Reining, V. Olevano, A. Rubio, and G. Onida, Phys. 
Rev. Lett. 88, 066404 (2002). 

[68] S. Botti, E. Sottile, N. Vast, V. Olevano, L. Reining, H.- 
C. Weissker, A. Rubio, G. Onida, R. Del Sole, and R. W. 
Godby, Phys. Rev. B 69, 155112 (2004). 

[69] A. Togo, E. Oba, and 1. Tanaka, Phys. rev. B 78, 134106 
(2008). 

[70] C. E. Macrae, P. R. Edgington, P. McCabe, E. Pidcock, 
G. P. Shields, R. Taylor, M. Towler, and J. van de Streek, 
Journal of Applied Crystallography 39, 453 (2006). 

[71] K. A. Udachin, C. 1. Ratcliffe, G. D. Enright, and J. A. 
Ripmeester, Supramolecular Chemistry 8, 173 (1997). 

[72] L. Pizzagalli, EPL 104, 56005 (2013). 

[73] S. Botti, M. Amsler, J. A. Elores-Livas, P. Ceria, 
S. Goedecker, and M. A. L. Marques, Phys. Rev. B 88, 
014102 (2013). 

[74] A. Bautista-Hernandez, T. Rangel, A. H. Romero, G.-M. 
Rignanese, M. Salazar-Villanueva, and E. Chigo-Anota, 
Journal of Applied Physics 113, 193504 (2013). 

[75] L. Genovese, A. Neelov, S. Goedecker, T. Deutsch, 
S. A. Ghasemi, A. Willand, D. Caliste, O. Zilberberg, 
M. Rayson, A. Bergman, and R. Schneider, J. Chem. 
Phys. 129, 014109 (2008). 

[76] C. Hartwigsen, S. Goedecker, and J. Hutter, Phys. Rev. 
B 58, 3641 (1998). 

[77] ASTM Standard G173-03(2008), Standard Tables for 
Reference Solar Spectral Irradiances: Direct Normal and 
Hemispherical on 37° Tilted Surface (ASTM Interna¬ 
tional, West Conshohocken, PA, 2008). 



9 


TABLE L Structural data of the allotropes containing cage-like building blocks. Column 1 contains our structural label, followed 
by the space group in the Hermann-Mauguin notation with superscripts indicating the number of atoms in the conventional 
cell. Column 3 contains the characterization of the building blocks, with the different maximum guest molecule radii in the 
cavities in column 4. Columns 5, 6, and, 7 contain the energy per atom with respect to d-Si, the volume per atom and the bulk 
modulus, respectively. The last column contains the corrected PBE band-gap. “D” indicates a direct or quasi-direct band-gap. 

Label Space-group Building blocks Guest radii (A) Energy (eV) Vol (A^) Bq Gap (eV) 
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FIG. 4. Energies as a function of volume/atom of the structures C01-C44 with respect to d-Si, together with the silicon 
clathrates type I and type II, and the predicted silicon allotrope Si20T which exhibits the best absorption to date [11]. 
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FIG. 5. Absorption spectra of selected phases. The reference air mass 1.5 solar spectral irradiance is shown in yellow (gray 
solid line), whereas the absorption of d-Si is shown by the solid black line. 
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